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ABSTRACT 



Aims. We present accurate photometric time series of two Be stars: NW Ser and V1446 Aql. Both stars were observed at the 
Observatorio de Sierra Nevada (Granada) in July 2003 with an automatic four-channel Stromgren photometer. We also present a 
preliminary theoretical study showing that the periodic variations exhibited by these stars can be due to pulsation. 
Methods. An exhaustive Fourier analysis together with a least-square fitting has been carried out on the time series for all four 
Stromgren bands. Several independent frequencies and non-periodic trends explain most of the variance. A theoretical non-adiabatic 
code applied to stellar models for these stars shows that g-modes are unstable. 

Results. Both stars show rapid variations in amplitude, probably due to a beating phenomenon. Four significant frequencies have 
been detected for each star. Comparison of the observed amplitude ratios for each pulsational frequency with those calculated from 
theoretical pulsation codes allows us to estimate the pulsation modes associated with the different detected frequencies. NW Ser seems 
also to show unstable p-modes and thus could be one of the newly discovered fi Cephei and SPB hybrid stars. Further spectroscopic 
observations are planned to study the stability of the detected frequencies. 

Key words, stars: oscillations (including pulsations) - stars: emission-line, Be - stars: individual: NW Ser, V1446 Aql 



1. Introduction 

Short-term variability, with periods ranging from a few hours to 
two days, is present in most of the early-type Be stars and in 
a significant fraction of late Be stars. It has been detected both 
spectroscopically, in the form of line profile variability (Ipv) and 
photometri cally, by means of non-regular light curves. 

iBaadd (I1982I) attributed the short periodic Ipv to non-radial 
pulsations (nrp). Further detailed studies of a few well-observed 
objects led to a complet e and consistent multimode modelling of 
the observed varia tions (Rivinius et al. 2001; Ma intz et al.l2 003; 
iNeiner eta l. 2005). In addition, high-precision photometric data 
obtained with the MOST satellite showed the presence of multi- 
ple perio ds in the Be stars ^-O ph dWalker et al.ll2005bl) and HD 
163 868 dWalker et al.ll2005al) . They modelled the detected fre- 
quencies as nrp in terms of p-modes in the case of f-Oph, and 
g- and r-modes in the case of HD 163 868. 

Similarly, COROT will provide important clues to under- 
stand the link between the beating of nrp modes and the origin 
of the Be phenomenon. COROT is based on ultra high preci- 
sion, wide field photometry for very long continu ous observing 
runs in the same field of view dBaglin et al.ll2.Q02b . In prepara- 
tion for this mission, a total of 84 Be stars have been studied 
photometrically in order to select the most suitable candidates to 
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be observed by COROT (Gutierrez-Soto et al., in preparation). 
Here we present the analysis of the photometric time series for 
two Be stars, namely NW Ser and V1446 Aql. 

NW Ser (HR 6873, HD 168797, V=6.14) is a bright and 
extensively observed B2.5IIIe star, which has been photomet- 
rically monitored for over two decades. From the analysis of the 
Hipparcos photometric data, severa l authors have obtained simi- 
lar results: iHubert & Floquej dl998l) found short- term variability 
with a period of 0.488 days; iPercy et all (Qj999) reanalysed the 
Hipparcos data together with ground-based photometry finding 
a similar period of 0.46 days, as well as a longer period of 5.5 
days; and f inally, two pe riods (0.475 and 0.406 days) have been 
detected by|Aerts|(2000) using only Hipparcos data, although the 
author claimed that this finding should be confirmed by means 
of ground-based observations. 

V1446 Aql (HD 179405, B2IVe) has been observed to be 
an emission -line star in the Mount W ilson objective prism all- 
sky survey (iMerrill & Burwelllll943l) . It appears as an irregular 
variable star fr om the Hipparcos data, which led to its inclusion 
in the GCVS dKazarovets et alj fl999). although no short-term 
periodic variability had been detected so far. 



2. Observations and frequency analysis 

Observations were made with the 0.9 m telescope of the 
Observatorio de Sierra Nevada (OSN, Granada, Spain) between 
July 1 and 9 in 2003. The instrument used is an automatic four- 
channel photometer, which allows simultaneous observations 
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Table 1. Comparison and check stars used in the differential pho- 
tometry. 





Star 


V 


Spectral Type 


variable 


NW Ser 


6.14 


B2.5IIIe 


comparison 


HD 170200 


5.71 


B8III-IV 


check 


SAO 123 607 


8.6 


B8 


variable 


V1446 Aql 


9.12 


B2IVe 


comparison 


HD 179 846 


8.29 


B8 


check 


HD 178 598 


9.45 


B8 



Table 2. Photometric precision in the instrumental system. 



Star 



u filter v filter b filter y filter 



NWSer 
V1446 Aql 



0.010 
0.011 



0.005 
0.007 



0.005 
0.007 



0.005 
0.008 



through the four uvby filters of the Stromgren photometric sys- 
tem. The data discussed here and presented in Tables A. 1 and 
A. 2 are the differential magnitude in the instrumental system be- 
tween the variable and the comparison star of Table Q] A check 
star has also been observed to verify possible intrinsic variations 
of the comparison star. Data have been corrected for sky back- 
ground and atmospheric extinction. Table[2]shows the photomet- 
ric precision as given by the standard deviation of the difference 
between the comparison and check stars for the whole campaign 
for each filter. The light curves of the stars are depicted in Figs.Q] 
and [3] only the v filter data are shown for clarity. 

The period analysis has been performed by means of stan- 
dard Fou rier analysis and le ast-square fitting. We have used 
Period04 (Lenz & Breger 2005) which is especially designed for 
the analysis of time series containing gaps. This program finds 
the frequencies one by one by computing the Fourier Transform 
and then adjusts the parameters of a sinusoidal function using a 
least-square fitting. This frequency is then removed and a new 
step is started finding a new frequency, the subsequent least- 
square fitting is performed allowing the two frequencies to move 
in order to obtain the minimum variance. The method is iterative 
and stops when removing of a new frequency is not statistically 
significant. 

We also used a non-linear multi-param eter fitting code which 
scans a wide range in frequency based o n I Vamcekl (1 1 97 lb and is 
explained in detail in lZerbi et al.l (1 1997b . This code is also well- 
suited for our case, for which daily aliases are present in the 
periodogram, due to the fact that our observations have been ob- 
tained at only one site. 

The way in which we determine whethe r the frequencies are 
statistically significant or not is described in Bre ger et~ai1 (I 1993b 
and basically consists of the calculation of the signal to noise 
ratio (SNR), the noise being the average amplitude, within a 5 
c d _1 frequency interval, of the residual periodogram after the 
prewhitening as previously explained around a frequency and 
the signal be i ng the corresponding amplitude of that frequency. 
iBreger et al.l (LL993) showed that this value must be greater than 
4. 

Frequencies, amplitudes and phases obtained for NW Ser 
and V1446 Aql are presented in Tables [3] and ^respectively for 
the four Stromgren bands. We also show the total fraction of the 
variance removed from the signal (R), the SNR and the <x of the 
final residual o- res . In Figs.|2]and|4]we display the successive pe- 
riodograms and the spectral window for NW Ser and V1446 Aql 
respectively, results are shown only for the v filter data for clar- 
ity. 
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Fig. 1. Light curve of NW Ser in the v filter. 



Following Montgomery & O'Donoghuel (0999), the ex- 
pected error in frequency for uncorrected observations can be 
derived from the equation 



V6 



CTf 



A ■ y/N ■ T 



where A/cr„ indicates the SNR, N is the number of obser- 
vations and T the time elapsed between the first an d the last 
data point. As noted by Schwarzenber g-Czernvl ( 1 1 99 lb . correla- 
tions in the residuals of the fitting have to be taken into account 
multiplying the error frequency by Vd, where D is the corre- 
lation length. D can be estimated by performing an autocorre- 
lation analysis of the final residuals. In our case the correlation 
length is 4. Therefore, we derive a formal error in frequency of 
6 - 20 x 10~ 4 c d l on frequencies Fl to F4 respectively for 
NW Ser and 8 - 16 x 10~ 4 c d -1 for V1446 Aql. 



2.1. NWSer 

The rapid amplitude variation over the short timescale spanned 
by our observations strongly suggests the presence of a beating 
phenomenon. This would explain why after prewhitening for the 
first peak (seen at Fl = 1.197 c d -1 in Fig. [2] together with its 
aliases at 1,2,... c d~') a second peak appears at a very close po- 
sition corresponding to a frequency at F2 = 1.126 c d -1 (together 
with its daily aliases). Note that the width of the sidelobe is al- 
most twice that calculated for the spectral window. This would 
mean that another frequency close to 1 . 197 c d -1 is also present. 
Following lLoumos & Deemingl (fT978). this occurs when the true 
frequencies differ by less than the half-width of the sidelobe and, 
in our case, a peak at the midpoint of the true frequencies would 
appear in the periodogram, as seen in Fig. |2] Fl and F2 have 
comparable amplitudes, which would produce the rapid ampli- 
tude variability in terms of a beating phenomenon. 

These two frequencies differ by less than the Rayleigh reso- 
lution, which is in our case 0.11 c d -1 . This could indicate that 
the value we calculate for these frequencies might not be the 
real ones. However, the rapid amplitude variation present in the 
light curve, could hardly be explained by any other origin than 
the presence of two closely spaced frequencies, even if their real 
values are not exactly those given in Table [3] Therefore we con- 
clude that two close frequencies do actually exist, although their 
exact values might differ f rom those given here. In any case, 
lLoumos & Deeming! dl978b noted that the true frequencies al- 
ways fall at the same midpoint as those found when the resolu- 
tion is not high enough. 
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Fig. 2. Successive periodograms for NW Ser for the v filter. The ticks indicate the positions of the detected frequencies and the 
dashed horizontal line indicates the Act level after the final prewhitening. 



Table 3. Frequencies, amplitudes and phases obtained with the 
multi-parameter fitting code for NW Ser for the four Stromgren 
bands. The SNR, the <x of the final residual (cr res ) and the total 
fraction of the variance removed from the signal (R) are also 
given. 



No. 


Freq. 


Amp. 


Phase 


SNR 


CTres 


R 




cd- 1 


mmag 


2;rrad 




mmag 


% 




Filter y 




N=162 




= 34.8 mmag 




Fl 


1.190 


38.1 


0.84 


25 


20.8 




F2 


1.119 


21.8 


0.88 


15 


12.1 




F3 


3.294 


13.6 


0.38 


7 


8.7 




F4 


1.415 


7.8 


0.29 


5 


7.1 


96 




Filter b 




N=162 




= 34.9 mmag 




Fl 


1.196 


36.6 


0.85 


25 


22.0 




F2 


1.127 


26.9 


0.86 


19 


11.7 




F3 


3.301 


11.6 


0.36 


7 


8.7 




F4 


1.414 


9.5 


0.28 


6 


6.2 


97 




Filter v 




N=162 


0"init 


= 35.1 mmag 




Fl 


1.197 


36.4 


0.86 


28 


22.6 




F2 


1.126 


27.7 


0.88 


22 


12.1 




F3 


3.304 


11.3 


0.30 


7 


8.5 




F4 


1.412 


10.5 


0.36 


7 


6.2 


97 




Filter u 




N=162 




= 44.4 mmag 




Fl 


1.159 


55.6 


0.95 


38 


26.2 




F2 


1.103 


33.3 


0.83 


23 


14.2 




F3 


2.274 


17.5 


0.49 


11 


10.5 




F4 


2.387 


12.9 


0.34 


8 


7.4 


97 



Residuals after prewhitening for these two frequencies have 
been analysed for additional periodic components. Frequencies 
F3 = 3.304 c d -1 and F4 = 1.412 c d _1 and their daily aliases 
have been detected in successive periodograms. These new fre- 
quencies have a much lower amplitude than Fl and F2, although 
they have a SNR greater than 4, as shown in Fig. [2] Note that 
F3 = 3.304 c d~ could be a combination of the two high- 
amplitude frequencies Fl and F2, since F3 ~ Fl +F2+ 1. Finally, 
the four frequencies which minimise the residuals are those 
shown in Table[3]for the four Stromgren bands. No other signif- 



icant peaks appear in the periodogram after the final prewhiten- 
ing. 

It is important to note that the four frequencies obtained in- 
dependently for the different vby filters are the same within the 
error boxes. For the u filter we found the led 1 alias of the lower 
amplitude frequencies F3 and F4. This is due to the low SNR of 
the data at that wavelength. As seen in Table[3] the percentage of 
the total fraction of the variance that is removed (R-values) for 
all filters is very high. 

Note that the final frequencies can be contaminated by 
1 c d -1 aliases, since the observations have been obtained at only 
one site. Keeping this in mind, if we compare with the frequen- 
cies found by other authors, the frequency F2 could be a 1-day 
alias o f the first frequency (f 1 = 2.11 c d~') detected by lAertsI 
(2000) using Hipparcos data and F4 could be a 1-day alias of 
her second frequency (f2 = 2.46 cd -1 ). In the same way, Fl 
co uld be the 1-day a lias of the frequency 2.17 cf ' obtained 
bv lPercv et al.l d!999h . The complexity of the Hipparcos spectral 
window due to the few (76) datapoints spanning ~ 1000 days 
makes the detection of multiple short periods very difficult and 
uncertain. Nevertheless, to test our results we performed a re- 
analysis of the Hipparcos light curve with the non-linear least- 
sq uare fi t ting c ode. We have detected the two frequencies found 
bv lAertsI (j2000) as the best 2-frequency fitting. However, the best 
3-frequency model is given by the frequencies 2.09, 1.43 and 
3.30 c d~'. The frequency F2 could be the 1-day alias of 2.09 
c d -1 and the two latter are similar to F4 and F3 respectively. 

In Table [3] we list the frequencies resulting from the anal- 
ysis of the OSN data. The above considerations indicate that 
these frequencies can be either real or in some cases their 1- 
day aliases. The analysis of the Hipparcos data shows the pres- 
ence of the same frequencies or their aliases, confirmin g that 
NW S er is a multiperiodic variable, as first suggested by Aerts 
(2000), and that the pulsation modes remain stable at least be- 
tween the two epochs. To ascertain what frequencies are real or 
the aliases, a more extended photometric dataset with a better 
spectral window, either from a ground-based multisite campaign 
or from space, would be needed. 
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Fig. 3. Light curve of V 1446 Aql in the v filter. 
Table 4. As for Tabled but for V1446 Aql. 



No. 


Freq. 


Amp. 


Phase 


SNR 


@~ res 


R 




cd- 1 


mmag 


2nmd 




mmag 


% 




Filter y 




N=136 


0" in it = 


23.89 mmag 




Fl 


1.631 


21.1 


0.76 


17 


19.6 




F2 


2.792 


18.8 


0.21 


11 


15.3 




F3 


2.559 


18.7 


0.66 


12 


12.3 




F4 


1.264 


13.0 


0.32 


11 


8.8 


87 




Filter b 




N=136 


@~ in it = 


24.85 mmag 




Fl 


1.622 


21.7 


0.78 


14 


20.4 




F2 


2.783 


20.9 


0.22 


11 


15.1 




F3 


2.555 


17.7 


0.65 


10 


12.3 




F4 


1.266 


12.9 


0.33 


8 


8.9 


87 




Filter v 




N=136 


Cmit = 


24.84 mmag 




Fl 


1.617 


22.1 


0.79 


15 


19.9 




F2 


2.783 


20.9 


0.23 


11 


14.3 




F3 


2.557 


17.5 


0.64 


10 


11.7 




F4 


1.269 


12.4 


0.32 


8 


8.3 


89 




Filter u 




N=136 


& init = 


29.10 mmag 




Fl 


1.601 


27.0 


0.77 


13 


22.6 




F2 


2.759 


25.5 


0.32 


11 


17.4 




F3 


2.514 


11.1 


0.67 


5 


14.6 




F4 


1.187 


15.5 


0.62 


8 


11.1 


85 



2.2. V 1446 Aql 

The light curve of V1446 Aql, presented in Fig. [3] also shows 
rapid variation which is probably caused by the presence of a 
beating phenomenon between two close frequencies. As the fre- 
quencies found for the four Stromgren filters are basically the 
same, we only discuss here the detailed analysis for the v filter, 
since it shows the highest SNR. Frequencies at 1.617 c d and 
at 2.783 c d~' and their daily aliases appear clearly in the pe- 
riodogram. Note that these two frequencies have similar ampli- 
tude. The amplitude of the first frequency (1.617 c d~') fulfils the 
signal to noise criterion, as shown in Fig. [4] After prewhitening 
for this frequency, the second frequency at 2.783 c d , and its 
daily aliases, appear clearly. By prewhitening for these two fre- 
quencies, other components have been found. Either a frequency 
at 2.557 c d -1 or its 1-day alias at 1.557 c d can be considered 
as the third frequency, although 1.557 c d -1 is so close to Fl that 
the fitting algorithm does not converge and we obtain unrealistic 
values for this solution. Therefore we remove the frequency at 
F3 = 2.557 and identify the fourth component F4 = 1 .269 c d -1 , 
or one of its aliases. 

The final set of frequencies, displayed in Table [4] fulfils the 
signal to noise requirement in all the filters. We stopped the fre- 
quency search at this point because new peaks did not fulfil the 
SNR criterion mentioned above. 



As in the previous subsection, we list the frequencies corre- 
sponding to the best fit to our OSN data. Due to the effects of 
data sampling discussed above, some of the listed frequencies 
might be 1-day aliases of the real ones. 



3. Theoretical modelling 

To determine whether the observed frequencies are predicted un- 
stable or not by the current models of stellar pulsation we have 
made a preliminary theoretical study. The physical parameters of 
NW Ser and V1446 Aql have been obtained from the analysis of 
high resolution spectra dFremat et al.ll20061) . taking into account 
the effects of gravitational darkening and assuming bo th stars are 
rotati ng at 88% of their critical break-up velocity (Frema t et al.l 
120051) . The obtained values and associated errors are depicted in 
Fig- El Note that both error boxes are located in th e overlapping 
regio n of the fi Cephei and SPB instability strips dPamvat nvkh 
119991) . 

Evolut ionary track s computed with the numerical code 
CESAM dMorellll997l) are shown in Fig. [5] First-order effects 
of rotation are taken into account in the equilibrium models. To 
do so, the equilibrium equatio ns are modified in the CESAM 
code in the manner described in Kippenhahn & Weigert ( 1990). 
The so-called pseMc/o-rotating models include the spherically av- 
eraged contribution of the centrifugal acceleration by means of 
an effective gravity geg = g — tft c (r), where g is the local gravity, 
r is the radius, and 3i c {r) — | rQ 2 (r) is the centrifugal acceler- 
ation of matter elements. This spherically averaged component 
of the centrifugal acceleration does not change the order of the 
hydrostatic equilibrium equations. The models are assumed to 
rotate uniformly and to conserve their total angular momentum 
during evolution. 

Considering the corresponding error boxes, masses in the 
range of 8.5-9.5 M are found to be representative of NW Ser. 
Rotational velocities are in the range of 250 to 270 km s -1 , radii 
from 5.8 to 7.9 R Q and ages around 20 Myr. In the same way, for 
V1446 Aql, a mass range of 7-7.5 M is found, with rotational 
velocities ranging from 238 to 33 1 km s , radii from 4 to 4.5 R 
and an age that brings the star close to the ZAMS. 

The in stability analysis i s then performed using the GraCo 
code (see iMova et aUl2004l for more deta ils), which i s base d 
on the non-adiabatic equations derived by lUnno et a 

D dl989l) . 

For both stars, the £ — 2 and i — 3 modes are predicted to be 
unstable. In the case of NW Ser, the SPB and /3 Cephei instability 
ranges are very close (Pig. [6j, and this suggests that NW Ser is a 
good candidate to be a hybrid pulsator. This result is consistent 
with its location in the HR diagram region where the instability 
strips of SPB and [} Cephei overlap. In the case of V1446 Aql, 
only high-order g-modes with 1-2 and 3 are predicted to be 
unstable close to the observed frequencies, as shown in Fig. [7] 

In order to have an idea of what spherical degree I is excited 
in these stars we ha ve appli e d the photometric mode identifica- 
tion as describ ed in I Watson] dl988l) and used for 6 Scuti stars in 
iGarridol (|2000|). We have computed the amplitude ratios for each 
pulsational frequency and we have compared them with those 
predicted for stellar models using different ^-degrees. In Figs. [8] 
and [9] we show the amplitude ratios for both stars. In the case 
of NW Ser, pulsation modes of degree I — 3 appear to be more 
probable for frequencies F2 and F4, whereas for Fl and F3 the 
more probable degrees seem to be i — 1 and 2. For V1446 Aql 
all detected frequencies seem to be associated with dipoles and 
quadrupoles. The u amplitude for frequency F3 is very low and 
probably due to the noise of the data at this wavelength. 
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Fig. 4. Successive periodograms for V1446 Aql for the v filter. Symbols are the same as for Fig. [2] 




4.30 



4.25 



4.20 



Fig. 5. Hertzsprung-Russel diagram showing the two Be stars 
considered in this study and some relevant evolutionary tracks 
selected for the modelling. The empty and shaded regions repre- 
sent the coiTesponding error boxes for NW Ser and V1446 Aql 
respectively. Continuous and dotted lines represent evolutionary 
tracks of 7, 8 and 9.5 M© for rotating models and non-rotating 
models respectively. 



At this point only spectroscopic observations can supply an 
unambiguous mode identification for these rapidly rotating ob- 
jects. 



4. Conclusions 

Accurate photometric light curves of the Be stars NW Ser 
and VI 446 Aql obtained during a preparatory program for the 
COROT space mission have been analysed. Our spectral anal- 
ysis reveals the presence of four statistically significant fre- 
quencies in each star. Together with the multiperi odicity of HP 
163 86 8 and £-Oph recently observed from space dWalker et al.l 
2005a, b), the present results point towards the interpretation of 
short-period variability of Be stars in terms of nrp, as suggested 
by spectroscopic observations obtained in the last decade. In this 




3 4 5 6 
Frequency (in cd" 1 ) 



Fig. 6. Growth rates diagram for NW Ser: positive values indi- 
cate unstable modes and viceversa. 
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it 
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Fig. 7. Growth rates diagram for VI 446 Aql. 



framework, a theoretical instability analysis of the observed fre- 
quencies has been perfo rmed. Recent studies dSaio et alj 2007; 
iDziembowski e t al. 2007) show that, for very rapid rotators, such 
as Be stars, the effect of rotation may also affect the instability 
of the modes. In future studies these effects should then be taken 
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Fig. 8. Observed amplitude ratios Ax/A y and their error boxes 
for the detected frequencies of NW Ser. Ax stands for any of 
the amplitudes in the uvby Stromgren filters. The lines represent 
the theoretical non-adiabatic predictions of the stellar models ex- 
plained in the text. 
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Fig. 9. Same as Fig.E but for V1446 Aql. 
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into account and compared with the present results. It is found 
that NW Ser could be an hybrid pulsator, since the ranges of un- 
stable g-modes (SPB-like) and p-modes (J3 Cephei-like) are very 
close. This would thus imply that NW Ser is a good candidate 
of the very few stars detected so far showing such pulsational 
characteristics. In the case of V1446 Aql, only frequencies cor- 
responding to predicted g-modes have been detected. 

Longer timebase spectroscopic and photometric observa- 
tions are thus required in order to confirm these results. 
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